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ABSTRACT 

We report on XMM-Newton spectroscopic observations of the luminous, 
radio-quiet quasar PDS 456. The hard X-ray spectrum of PDS 456 shows a 
deep absorption trough (constituting 50% of the continuum) at energies above 7 
keV in the quasar rest frame, which can be attributed to a series of blue-shifted 
K-shell absorption edges due to highly ionized iron. The higher resolution soft 
X-ray RGS spectrum exhibits a broad absorption line feature near 1 keV, which 
can be modeled by a blend of L-shell transitions from highly ionized iron (Fe xvil 
- xxiv). An extreme outflow velocity of ~ 50000 km s~^ is required to model 
the K and L shell iron absorption present in the XMM-Newton data. Overall, a 
large column density {Nh = 5 x lO^^cm"^) of highly ionized gas (log ^ = 2.5) is 
required in PDS 456. A high mass outflow rate of ~ 10 Mq year~^ (assuming a 
conservative outflow covering factor of 0.1 steradian) is derived, which is of the 
same order as the overall mass accretion rate in PDS 456. The kinetic energy 
of the outflow represents a substantial fraction (~ 10%) of the quasar energy 
budget, whilst the large column and outflow velocity place PDS 456 towards the 
extreme end of the broad absorption line quasar population. 

Subject headings: galaxies: active — quasars: individual: PDS 456 — X-rays: 
galaxies 



1. Introduction 

PDS 456 is a luminous, low redshift {z = 0.184) radio-quiet quasar identified in 1997 
(Torres et al. 1997). The optical and infra-red spectra (Simpson et al. 1999) show broad 
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Balmer and Paschen lines (e.g. H/9 FWHM 3000 km s~^), strong Fe ll, a hard (de-reddened) 
optical continuum (/j, cx and one of the strongest 'big blue bumps' of any AGN 

(Simpson et al. 1999; Reeves et al. 2000). PDS 456 has a de-reddened, absolute blue mag- 
nitude of Mb ~ —27 (Simpson et al. 1999), making it at least as luminous as the radio-loud 
quasar 3C 273. Indeed PDS 456 is the most luminous known AGN in the local Universe 
(z < 0.3), its luminosity being more typical of quasars at z=2-3, at the peak of the quasar 
luminosity function. 

PDS 456 was first detected as the X-ray source RXS J172819.3-141600 in the ROSAT AW 
Sky Survey (Voges et al. 1999). Subsequent ASCA and RXTE and BeppoSAX ohscrvations 
showed that it was highly X-ray variable (Reeves et al. 2000, 2002), requiring that the X- 
ray source was very compact (< 3 Schwarzschild radii) and accreting with a unusually high 
efficiency (between 6% and 41%), close to the hmit for rotating Kerr black hole (Thorne 
1974). The X-ray spectrum of PDS 456, obtained by ASCA and RXTE also appeared to 
be unusual for a quasar. A strong edge like feature was present in the iron K-shell band, 
either due to a highly ionized absorber or reflector, whilst absorption towards the soft X-ray 
band was also observed (Reeves et al. 2000). An XMM-Newton observation of PDS 456 was 
subsequently obtained in AO-1, in order to obtain high signal to noise spectra of both the 
iron K-shell complex with the EPIC CCD detectors and a high resolution spectrum of the 
soft X-ray absorber with the XMM-Newton RGS (Reflection Grating Spectrometer). Here 
we report on the X-ray spectra obtained through the XMM-Newton observations of PDS 456, 
which reveal evidence for a massive outflow, with a velocity of ~ 50000 km s~^. A previous 
paper (Reeves et al. 2002) discussed in detail the X-ray variabihty during the XMM-Newton 
observation and a simultaneous week long BeppoSAX observation. 



2. The X-ray Observations 

PDS 456 was observed by XMM-Newton on 2001 February 26th with a duration of 40 
ks. Data were taken with the EPIC-pn detector (Struder et al. 2001) in Full Window Mode 
and with the EPIC-MOS detector (Turner et al. 2001) in Large Window mode. The data 
was reduced using version 5.4 of the XMM-SAS software. Note that the timc-avcraged, 0.5- 
10 keV flux of PDS 456 during the observation was 9.7 x 10~^^ ergs cm~^ s^^; at this flux 
level the effect of photon pile- up is small (< 1%). All time intervals were included, as the 
observation contained no background flares. Data were selected using event patterns 0-12 
(for the MOS) and pattern 0-4 (for the pn) and only good X-ray events were included. The 
spectra were extracted from circular source regions of 40" radius, whilst background spectra 
were extracted from an offset circle, close to PDS 456, but free of any background sources. 
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Response matrices and ancillary response files were generated using the SAS tasks rmfgen 
and ARFGEN respectively. 

Background subtracted spectra were fitted using XSPEC vll.2, including data over the 
energy ranges 0.4 to 12 keV for the EPIC-pn and 0.4 keV to 10 keV for the co-added EPIC- 
MOS spectrum. A Galactic absorption column of Nh — 2x 10^^ cm~^ (Dickey & Lockman 

1990) was included in all the fits. A single absorbed power-law fit (photon index, F = 2.0) 
to both MOS and FN data is clearly inadequate (reduced = 10). Figure 1 shows the 
data/model ratio to this power- law, a strong soft excess is observed below 1 keV and a large 
drop in counts is observed above 7 kcV in the iron K-shcll band. The MOS and pn spectra 
are consistent, except below 0.5 keV where there is a small divergence in spectral slope. 

3. The Ionized Absorber in PDS 456 

Initially, to model the iron K-shell band, the data were fitted from 2-12 kcV (or 2-10 kcV 
for EFIC-MOS), so that the effects of the large soft excess are negligible. As the EFlC-pn 
data contained significantly more counts than the EFIC-MOS at the highest energies, we 
proceeded to fit the EFIC-pn data only. The ratio spectrum above 2 keV to a power-law 
fitted over the 2-5 keV range (i.e. outside of the iron K band) is shown in the inset to Figure 
1. A strong deficit of counts occurs above 7 keV, where one would expect to observe K-shell 
absorption edges from highly ionized iron. Indeed a simple Galactic absorbed power-law 
produced a poor fit {x^/dof — 1001/746, null hypothesis probability 9 x 10~^°). 

Adding three absorption edges between 7 and 10 keV significantly improves the fit 
{x^/dof — 739/740), modehng the broad Fe absorption feature present in the X-ray spec- 
trum. The energies (in the quasar rest frame) and optical depths of the edges are then; 
^ = 9.3 ± 0.1 keV, E = 8.5± 0.2 kcV and E = 7.6 ± 0.3 kcV with r = 0.56 ± 0.11, 
r = 0.26 ±0.11, T = 0.24 ±0.08 respectively. The breadth of the absorption feature suggests 
that a wide range of ionizations may be present, whilst the depth of the feature implies a 
column of gas of Nh ~ 10^^ cm~^. Note that the upper limit to a narrow (i.e. unresolved) 
Fe emission line (in the energy range 6.4 - 6.97 keV) is < 12 eV, implying that the solid 
angle subtended by the absorber must be small {Q < 0.5 steradian). Note that an extreme 
broad iron Ka emission line from an accretion disk cannot model the shape of the spectrum 
from 2-12 keV. If one attempts to model the spectrum with a Fe fine profile from a rapidly 
rotating Kerr black hole (Laor 1991), the physical parameters required are unphysical (disk 
emissivity R^^^, inclination=90°, equivalent width > 1 kcV), because the line is (unsuccess- 
fully) attempting to model the broad absorption dip above 7 keV with the blue-wing of the 
line profile. 
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To obtain a more physical representation of the Fe K-shell absorber, we used a grid of 
models generated by the XSTAR photoionization code (Kallman et al. 1996) to fit the EPIC- 
pn data above 2 keV. Solar elemental abundances were assumed, whilst a turbulence velocity 
of 1000 km s~^ was used. A one zone photoionization model was adopted, with the column 
density (A^h), ionization parameter (^), partial covering fraction (/) and outflow velocity 
(^out) as free parameters. An effective hydrogen column density (for a solar abundance 
H/Fe ratio) of S.Tllg x 10^^ cm~^ is required, together with an ionization parameter of 
log ^ = 2.5 ± 0.3 and a partial covering fraction of / ~ 0.6, whilst the power-law continuum 
steepens to F = 2.40 ±0.05. At this ionization state most of the absorption is due to Fe xvil 

- Fe XXIV (i.e. where the L-shell is partially filled). The best fit spectrum modeled by 
XSTAR is shown in Figure 2 (panel a), the depth and energy of the absorption feature being 
well matched, whilst the partial covering warm absorber also models the convex spectral 
curvature (through the increase in opacity at lower energies due to iron L-shell absorption) 
present in the EPlC-pn data between 2-6 keV. Overall, the fit statistic for this model is 
excellent {x^/dof = 717/740). 

Importantly a large outflow velocity of 47000 km s~^ (90% confidence limits, 36000 km s~^ 

- 82000 km s~^) is needed to model the spectrum. Although the exact outflow velocity is not 
well constrained by the EPIC-pn data alone, a significant outfiow velocity is still required. 
Figure 2(b) illustrates this point, if one fixes the absorber velocity to zero, then the model 
fit becomes poor as the iron K edges then lie at a significantly lower energy than is observed 
in the EPIC data. Indeed if one tries to increase the energy of the model absorption feature 
by increasing the ionization state of the absorber (so that He and H-like iron dominates), 
the result is a worse fit still, partly because the absorption edge feature becomes weaker (as 
a significant fraction of iron becomes fully ionized) and also because the most prominent 
absorption features are now the strong Ka and Kf3 absorption lines due to Fe XXVI which 
are not present in the actual XMM-Newton data (Figure 2c). Thus a significant outflow 
velocity is required by the data, indeed one can place a conservative lower-limit (at 99.9% 
confidence) on the velocity of 24000 km (or ~ 0.08c). 



3.1. The soft X-ray RGS Spectrum of PDS 456 

Extrapolation of the hard X-ray spectrum to lower energies (down to 0.4 keV), leaves 
clear residuals in the EPIC MOS and pn data (e.g. Figure 1). A strong soft excess is 
present below 1 keV, whilst residuals due to an additional absorber may also be present in 
the soft X-ray band. To investigate this at higher spectral resolution, we used data available 
from the Refiection Grating Spectrometer (RGS) below 2 keV. The most noticeable feature 
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present is a broad absorption trough observed between 12-15 A (see Figure 3). The feature is 
highly significant, being detected at > 99.99% confidence. When parameterized by a simple 
broad Gaussian shaped absorption line, the best fit parameters (in the quasar rest frame) 
are E = 1.07 ± 0.02 keV with an equivalent width of 26 ± 10 eV and a velocity width of 
9000 ±3000 km s'^ 

As in the previous section, we used our grid of photoionization models to characterize 
the properties of the absorbing gas, apparent in the RGS spectrum. A lower column density 
(Ah ~ 5 X 10^^ cm^^) is required to fit the data, whilst both a similar ionization parameter 
(log ^ = 2.6 ± 0.5) and outfiow velocity (with ?;out = 570001^qqqq km s~^) are needed, as per 
the fits to the EPlC-pn data. Partial covering is not required in the RGS fits. The photon 
index obtained (F = 2.4, after correcting for absorption) is consistent with the EPIC-pn 
data, whilst a soft blackbody component (with kTsB = 80 ± 6 eV) is also needed in the fit. 
The lower column density derived from the RGS fit indicates that the absorber is patchy and 
thus is likely to be located quite close to the X-ray source. Thus we see a higher (partially 
covered) column along the line of sight towards the hard X-ray continuum, whilst the column 
towards the soft X-ray continuum is lower. 

Inspection of the best-fit photoionization model implies that a blend of L-shell transi- 
tions from Fe XVII - Fe XXIV as well as K-shell absorption from Ne ix/x is responsible for 
the broad absorption feature at 12-15 A. A weak Lj-a absorption line due to O VIII may 
also be present at 19 A. The outfiow velocity derived appears robust, due to the positioning 
of the broad absorption feature. If one fixes the absorber velocity at zero in the quasar rest- 
frame, then the Fe L-shell absorption in the model is located at energies below the observed 
absorption trough, whilst increasing the ionization of the absorber only serves to significantly 
weaken the soft X-ray absorber. 



4. Discussion 

X-ray spectroscopy of PDS 456 with XMM-Newton has revealed a high column density 
(Ah = 5 X 10^^ cm~^), highly ionized (log ^ = 2.5) warm absorber, its main signature being a 
series of K and L shell absorption edges and lines due to iron in the ionization range Fe xvil - 
Fe XXIV. Crucially, a near relativistic outflow velocity is required, of the order 50000 km s~^. 
This velocity is similar to the actual systemic velocity of PDS 456 (at z — 0.184), therefore it 
is possible that the absorbing material is unrelated to the quasar, and is local to our Galaxy. 
However we consider it extremely unlikely that such a high column, highly ionized absorber 
could occur by chance along the local line of sight. Furthermore, although the Galactic 
column towards PDS 456 is moderately high (A^h = 2 x lO^^cm"^), this neutral absorption 
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does not contribute towards the features seen in the iron K and L-shell bands, and in any 
event has been accounted for in all our modeling. Therefore it is highly probable that the 
absorber is intrinsic to the quasar PDS 456. 

4.1. A Massive Outflow in PDS 456 

The combination of a large column density, high ionization and a fast outflow velocity in 
PDS 456 indicates that the absorber is located close to the quasar's central engine and that 
its mass outflow rate is significant relative the to the accretion rate. A maximum distance 
to the absorber can be calculated on the condition that its thickness (AT?) cannot exceed 
its distance from the ionizing source (R), i.e. AR/R < 1. This requires that R < W^Rg 
(where Rg = GM/c^) for a black hole mass of lO^M© for PDS 456 (Reeves et al. 2000), 
thus it is located within the BLR. One likely possibility is that the absorber arises through a 
disk driven outflow or wind. The exact origin of these winds and the mechanism responsible 
for driving them (radiation and/or magnetic) is unclear. For Broad Absorption Line (BAL) 
QSOs, hue driven winds have been proposed (Murray et al. 1995; Proga, Stone & Kallman 
2000) as well as centrifugal driving along open magnetic field lines (Bottorff et al. 1997). 

Using conservation of mass, the outflow rate (Mout) for an outflow subtending a solid 
angle of Q steradian with a constant velocity I'out is given by:- 

Mout = VtnPvontm^ 

where n is the density, r the distance from the X-ray source and mp the proton mass. 
Now the ionization parameter ^ is defined as ^ = Lx/nr^, so we can obtain:- 

Mout = QvoutrripLx/^ 

Prom the observations, the hard X-ray luminosity of PDS 456 (above 2 keV) is ~ 
10^^ erg s~^, the ionization parameter is ^ ~ 1000, whilst v — 50000 km s~^. Prom the lack of 
any iron K line emission it is likely that the outfiow subtends a solid angle of < 0.5 steradian. 

However even for a conservative value of 0.1 steradian, the mass outflow rate derived is 
Mout = 8 X 10^^ g s~^ or ~ lOM© year~^. A lower limit of ~ 5Mq year"-*^ is obtained using 
the 99.9% confidence lower bound on the outfiow velocity of Vont = 24000 km s^^ Note the 
outfiow rate calculated assumes the gas is in equlibrium, if the gas is being accelerated then 
the actual rate may differ. 

PDS 456 is an extremely luminous quasar with a bolometric luminosity of L^oi ~ 
]^q47 gj,g g-i ^^^^ ^ black hole mass of IO^Mq assuming accretion occurs near the Eddington 
limit (Simpson et al. 1999; Reeves et al. 2000). The mass accretion rate for PDS 456, given 
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hy M — Lboi/rjc^, is ~ 20Mq year""^, for an efficiency oi rj — 0.1. Tfius it seems likely the 
outflow rate is of the same order as the mass accretion rate in PDS 456. Similarly the kinetic 
output of the outflow in PDS 456 is large, about 10^^ erg s~^ corresponding to 10% of the 
bolometric output of the quasar. The high velocity also implies that the outflow is likely to 
be driven from the accretion disk at a fairly small radius. For a Keplerian disk, the launch 
radius for the outflow is simply i?out = hlsc^gi where v^sc is the escape velocity of the 
material. Assuming that v^sc ~ ^out (although some acceleration of the matter will probably 
occur), this imphes that the outflow in PDS 456 is launched from a radius of ~ 40i?g (about 
10^^ cm for a lO^M© black hole). Indeed the launch radius may be smaller if the gas is not 
in simple equilibrium. 

PDS 456 appears to be in a small but growing group of AGN to show extreme absorption 
properties in the iron K-shell band. The lensed BAL quasar APM 08279+5255 has been 
found to exhibit iron K absorption hues and/or edges (Chartas et al. 2002; Hasinger, Schartel 
& Komossa 2002) in both Chandra and XMM-Newton observations. A similar absorber has 
also been reported in the BAL quasar PG 1115+080 (Ghartas. Brandt & Gallagher 2003). 
The outflow velocities are similar to PDS 456, in the range 0.1c — 0.4c. Given this apparent 
association between high velocity X-ray outflowing gas and the UV BALs, one may predict 
to see BAL features in the UV spectrum of PDS 456. Indeed a HST-STIS spectrum of PDS 
456 (O'Brien et al. 2003, in preparation) reveals that there may be cooler gas in the outflow; 
a BAL is seen in Lyman-cc with an outflow velocity of 12000 - 22000 km s~^, whilst the G 
IV emission line is also blue-shifted by ~ 5000 km s~^. Thus PDS 456 may just reside at the 
extreme end of the BAL quasar population (Turnshek et al. 1988; Weymann et al. 1991). 

Similar highly ionized X-ray outflows have recently been discovered in the non-BAL 
radio-quiet quasars PG 1211+143 and PG 0844+349 (Pounds et al. 2003a,b), suggesting 
that these properties are not unique to UV identifled BAL quasars. We also note that 
strong iron K absorption has also been observed in the NLSls IH 0707-495 and IRAS 13324- 
3809 (BoUcr ct al. 2002, 2003) and also in some X-ray binaries, e.g. GRS 1915+105 (Lee 
et al. 2002), GX 13+1 (Sidoli ct al. 2002), X 1624-490 (Parmar et al. 2002). In the quasars 
PDS 456, PG 1211+143 and PG 0844+349, the X-ray columns are near to lO^^cm^^ the 
outflow velocities measured are ~ 0.1c, whilst the ionizations are also similar ~ 10^). 
Using PG 1211+143 as an example. King & Pounds (2003) have recently demonstrated, for 
radiatively driven outflows, that objects accreting near to the Eddington rate are likely to 
contain massive (Mout ~ Medd) Gompton-thick winds. Indeed when the outflow rate is of 
the same order as the accretion rate, the flow can become optically thick at ~ 10 — lOORg. 
This can account for both the high columns and velocities of the X-ray outflows, as well as 
provide a natural source of thermalized photons for the UV and soft X-ray excesses observed 
in these quasars. 
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However other scenarios are possible, for instance Reeves et al. (2002) recently used a 
magnetic flare model to explain the rapid X-ray variability in PDS 456, through coherent 
flaring events occuring when the accretion rate is close to Eddington. Indeed it is conceivable 
that the outflowing matter in PDS 456 is the ejecta associated with these flares. Interestingly 
the output of the magnetic events (~ 10'^^ erg s~^) in PDS 456 is very similar to the kinetic 
power of the outflow, so this scenario seems energetically plausible. Overall the extreme 
properties of the outflow in PDS 456 seem hnked to its high accretion rate (Reeves et al. 
2000), which appears to drive the energetics of its unusually massive outflow. 
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Fig. 1. — Ratio of the broad band X-ray spectrum of PDS 456 to a simple Galactic absorbed 
power-law model (with T = 2.0). EPIC-pn data are plotted in black, EPIC-MOS in gray. A 
clear soft X-ray excess is observed below 1 keV, as well as a large decrease in counts above 
7 keV. The inset plots the data/model ratio to a power-law fit to the 2-12 keV EPIC-pn 
spectrum in greater detail. A 50% drop in counts above the iron K-shell band is observed. 
Note the energy scale is plotted in keV in the observed frame. 
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Fig. 2. — Warm absorber fits to the EPIC-pn data of PDS 456 using tlie xstar pliotoioniza- 
tion code, (a) sliows the best fit model (solid line) fitted to the EPIC-pn data (crosses), which 
requires an outflow velocity of ~ 50000 km s~^ and an ionization parameter of log ^ = 2.5; 

(b) shows the same model, but with no outflow velocity, resulting in a significantly worse fit; 

(c) the same model, but with no outflow and a higher ionization parameter of log ^ = 3.5. 
At this high ionization the model cannot reproduce the shape of the iron K-shell absorption 
seen in PDS 456. 
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Fig. 3. — XMM-Newton RGS-2 spectrum of PDS 456, the data points are shown as crosses. A 
strong absorption trough is present at 12-15 A, observed frame. The best fit photoionization 
model (sohd fine) is superimposed on top of the data points, a high outfiow velocity of 
~ 50000 km s^^ is required for the absorber. Note due to a malfunctioning chip, RGS-1 does 
not containing data in the 12-15 A range and thus is not plotted. Most of the opacity of the 
absorber results from L-shell transistions of iron in the ionization range Fe xvil - Fe xxiv. 



